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Combined results are presented on the search for a neutral Higgs boson in the di-tau final state 
using 1.8 fb -1 and 2.2 fb -1 of integrated luminosity collected at the CDF and DO experiments 
respectively. Data were collected in pp collisions at a centre of mass energy of 1.96 TeV during 
Runll of the Tevatron. Limits are set on the cross section x branching ratio ranging from 13.6 pb 
to 0.653 pb for Higgs masses from 90 GeV to 200 GeV respectively. The results are then interpreted 
as limits in four different benchmark scenarios within the framework of the MSSM. 



Preliminary Results 



* The Tevatron New Phenomena and Higgs working group can be contacted at TEVNPHWG@fnal.gov. More information can be found 
at |http: //te vnphwg.f nal. gov/| 
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I. INTRODUCTION 

Spontaneous symmetry breaking in the electroweak sector is an attractive solution to the problem of the origin 
of particle masses within the Standard Model (SM). However, extreme fine tuning is required to avoid divergencies 
in radiative corrections to the Higgs mass. Supersymmetry (SUSY) as an extension to the SM, provides a natural 
means to avoid this as well as potentially providinga candidate for dark matter and GUT-scale unification. The 
Minimal Supersymmetric Standard Model (MSSM) [l[ requires the introduction of two Higgs doublets and predicts 
the existence of five physical Higgs bosons after symmetry breaking: three neutral (h, H, and A) and two charged 
(H^). The ratio of the vacuum expectation values of the two doublets is denoted by tan/3. For high values of tan/3 two 
of the three neutral Higgs bosons have approximately the same mass and couplings. These couplings are enhanced 
with respect to the charged leptons and down-type quarks by a factor tan/3 relative to the SM, and suppressed for the 
neutrinos and up-type quarks. The near degeneracy contributes an additional factor two enhancement in the cross 
section. Thus for low Ma and high tan/3 the Tevatron can probe a number of benchmark scenarios in the MSSM 
complementing the regions of the SUSY parameter space probed by the LEP experiments 0- 

The results presented here represent an update to the previous combination Q. The same inputs, with an additional 
mass hypothesis at 90 GeV, are used but the treatment of correlations between the systematic uncertainties on those 
backgrounds estimated from Monte-Carlo simulations has been improved and the intepretation of the results in the 
MSSM uses calculations from the latest version of feynhiggsQ]. 

II. ANALYSIS SUMMARY 

The CDF and DO detectors are described in detail elsewhere [H, Q . The searches combined here are described in 
detail in 043 an d earlier published results from CDF and DO can be found in E[ ■ Searches are performed at CDF 
and DO for MSSM Higgs boson production with subsequent decays to taus in a number of channels characterised by 
the decay products of the r leptons. Included in this combination are 1.8 fb" 1 of data collected at CDF in three final 
states: r e Th a d, T^had and T e 7) t , (where r e , t m , and rtad denote r decays to electron, muon and hadrons respectively) and 
1.0 fb -1 in the same three channels and an additional 1.2 fb -1 in the r^rhad final state collected at DO. Additionally, 
the searches from DO are split further depending on the hadronic decay multiplicity. 

A. Lepton Identification 

Electrons are identified through their characteristic energy deposits in the calorimeters. Reconstructed clusters of 
energy in the calorimeters are required to be isolated and match a reconstructed track, suppressing photon back- 
grounds. Muons are identified by matching charged tracks in the central tracking detectors with hits in the muon 
detectors. Muon candidates are also required to be isolated in both the central tracking detectors and in the calorime- 
try. 

Hadronic decays of r leptons are identified at CDF by selecting isolated narrow clusters in the calorimeter with 1 
or 3 spatially matched charged tracks. These are reconstructed using a variable sized cone algorithm whose angle, a, 
is set to be the minimum of 10° and (5 GeV) / E c i radians, where E c i is the calorimeter cluster energy. Strict isolation 
limits on the number of tracks and the calorimeter energy within an annulus around the candidate from a out to an 
angle of 30° are used to suppress quark and gluon jets. In the case of three-prong candidates the sum of the charges 
of the tracks is required to be ±1. One-prong candidates are rejected if found to be consistent with an electron having 
undergone significant bremsstrahlung. 

In the DO analyses, the hadronic decays of the r are divided into three categories: r types 1 and 2 are one-prong 
candidates with energy either in only the hadron calorimeter (tt^ like) or in both the electromagnetic and hadron 
calorimeters (p^ like) respectively; r type 3 is a three-prong candidate with an invariant mass (constructed from the 
three tracks) below 1.7 GeV and matching energy deposits in the calorimeters. A neural network (NN) is trained 
for each type to separate hadronic tau decays from jets using MC Z — > tt as the signal and multi-jet events taken 
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from data as the background. An additional NN is trained on electron Monte-Carlo events and is employed to reduce 
backgrounds from electrons faking type 2 taus. 



The acceptance for signal is determined from Monte-Carlo simulations, using the pythia[12| event generator with 



bb — > A are considered by CDF, whereas at DO only gg — >• A is simulated - the acceptances are seen to be very similar 
for both production modes. In the interpretation of the results in the framework of the MSSM as limits in the tan/3- 
Ma plane both production modes are taken into account as well as an additional factor of approximately two on the 
cross section due to the near degeneracy of two of the three neutral Higgs bosons. Most Standard Model backgrounds 
have been generated with pythia: Z/7* — > , W -> tv, di-boson production, tt (comphep + PYTHIA) [16l| , W 
and Z boson samples where there is one or more additional jets in the final state have been simulated with ALPGEn[17| 
with matching to pythia for hadronization. Di-boson and tt samples are normalised using calculations to next-to-LO 
(NLO)m[ljfwhile Z/7* samples are generally normalised to next-to-NLO (NNLO) 

Events are selected by the trigger using inclusive electron and muon (DO) and lepton plus track (CDF) trig- 
gers and after offline reconstruction candidate events must contain two isolated opposite charged final state leptons 
(e, fj,, t). Leading sources of background are: — > tt, multi-jet, W — > ev, p,v,TV, Z — > fifi,Z — >ee, di-boson 
(WW,WZ, ZZ,Wj, Zj), and tt-pair production. In the t^t^ channel at CDF, events from the sidebands of the 
lepton isolation are used to determine the jet backgrounds. For the r e Th a d and T M Th a d channels the jet backgrounds 
where a jet fakes a Thud are estimated by weighting data events passing very loose cuts with the jet-r fake probability 
measured in an independent jet sample. The multi-jet contribution from data collected at DO is estimated using ei- 
ther T e Th a d candidate events where the electron and r have the same charge or using inverted lepton selection criteria 
(■7/jTiad and r e r M channels). The normalisation of the W production backgrounds is estimated from a data sample 
dominated by IV+jet events. 

In the TeThad and TeThad channels the electron or muon are required to be isolated and have a transverse momentum, 
Pt > 10 (CDF) or 15 (DO) GeV. One-prong hadronic tau candidates are accepted with p T > 15 GeV (CDF), 16.5 GeV 
(DO) and three-prong are required to have px > 20 GeV (CDF) 22 GeV (DO). Additional cuts are placed on the scalar 

sum of transverse momenta in the event at CDF, Ht = \p^ M | + |pTTh a d| + \ $t\ > 55 GeV, where $ T is the momentum 
imbalance in the transverse plane. In one-prong events where the rate at which jets fake taus is lower a slightly looser 
cut is used, Ht > 50 or 45 GeV for T e Th a d and TeThad respectively. Further cuts on the relative directions of the taus 



and the ^ T (CDF and DO) and the transverse mass (DO) Mt = y 2p T ,$ T (l — cosA</>), where Acj> is the azimuthal 
angle between the electron or muon and the hadronic tau, serve to suppress background contributions from VF+jets 
production. 

In the CDF t £ t^ channel events are selected requiring one central electron and one central muon with: mm(E T ,p T ) > 
6 GeV,max(E T ,p T ) > 10 GeV, and \E T \ + \p T \ > 30 GeV. DO make a similar selection, where: p T > 10 GeV and 
p T > 12 GeV and the invariant mass of the electron-muon pair exceeds 20 GeV and \E T \ + \p T \ + |$ T | > 65 GeV. Table 
U shows the expected number of backgrounds, observed events in data and the signal efficiency for Ma = 130 GeV. 

In setting the limits, events from regions of phase space with a similar ratio of expected signal (S) to background (B) 
can be combined without loss of sensitivity. Thus a useful way to visualize the comparison of expected backgrounds 
and the observed data is to show the event distributions binned in this ratio S/B. For the channels combined in 
the results presented in this note these distributions are shown in Figure [TJ The left hand plot is for a signal, 
Ma = 100 GeV and axBr= 2.0 pb and the right hand plot for a signal of Ma — 180 GeV and a x Br=0.66 pb. Good 
agreement is observed between the data and expected backgrounds. The integrals of these distributions starting from 
the high S/B side and working downwards are shown in Figure [21 displaying the signal +background, background-only 
and data sums. 



B. Signal, Backgrounds and Event Selection 
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CDF 






D0 




Source 


T e T M 


T e T had 






TeThad 


T M T had 


Z 


605 ± 51 


1378 ± 117 


1353 ± 116 


212 ± 20 


581 ± 5 


2153 ± 156 


Z ->■ e+e~/Z 


19.4 ± 5.7 


70 ± 10 


107 ± 13 


10.4 ± 1.3 


31 ± 2 


66 ± 8 


diboson + tt 


20.5 ± 7.0 


8.2 ± 4.2 


6.6 ± 3.7 


6.1 ± 0.6 


3.1 ± 0.3 


16 ± 3 


multi-jet + W->lv 


57.1 ± 13.5 


467 ± 73 


285 ± 46 


37.9 ± 7.7 


374 ± 48 


216 ± 41 


Total Background 


702 ± 55 


1922 ± 141 


1752 ± 129 


266 ± 22 


989 ± 82 


2451 ± 162 


Data 


726 


1979 


1666 


274 


1034 


2340 


Signal Efficiency /% 


0.32 ± 0.01 


0.77 ± 0.01 


0.67 ± 0.01 


0.41 ± 0.03 


0.73 ± 0.03 


0.99 ± 0.05 



TABLE I: Expected numbers of background and observed data events and signal efficiency for Ma = 130 GeV. Errors include 
full systematic uncertainties, that are in some cases correlated. 
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FIG. 1: Events binned by the ratio of expected signal to expected background for a signal of Ma = 100 GeV, and crxBr = 2.0pb 
(left) and M A = 180 GeV, and a x Br = 0.66pb (right) . 



III. COMBINATION 



To gain confidence that the final result does not depend on the details of the statistical formulation, two types 
of combinations are performed, using the Bayesian and Modified Frequentist approaches, which give similar results 
(within 10%). Both methods rely on distributions in the final discriminants, and not just on their single integrated 
values. Systematic uncertainties enter as uncertainties on the expected number of signal and background events, as 
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FIG. 2: Integrated distributions of S/B, starting at the high S/B side for M A = 100 GeV, crxBr= 2.0pb (left) and M A = 
180 GeV, a x Br=0.66pb (right). The total signal+background and background-only integrals are shown separately with data 
superimposed. Data points are only plotted for those bins with data events. 



well as on the distribution of the discriminants in each analysis ( "shape uncertainties" ) . Both methods use likelihood 
calculations based on Poisson probabilities. In all channels the visible mass distribution is used to set limits. 



A. Bayesian Method 



Because there is no experimental information on the production cross section for the Higgs boson, in the Bayesian 
technique [2l| a flat prior is assigned for the total number of selected Higgs events. For a given Higgs boson mass, the 
combined likelihood is a product of likelihoods for the individual channels, each of which is a product over histogram 
bins: 

Nc Nbins n n P 

C(R, s, b\n, 9) x tt(0) - I] II *Ce-"«/ny! x [] e"^ /2 (1) 

i=l j = l k=l 

where the first product is over the number of channels (Nc), and the second product is over histogram bins containing 
Uij events, binned in ranges of the final discriminants used for individual analyses, such as the di-jet mass, neural- 
network outputs, or matrix-element likelihoods. The parameters that contribute to the expected bin contents are 
fiij = R X Sij(6) + bij(6) for the channel i and the histogram bin j, where sy and bij represent the expected 
background and signal in the bin, and R is a scaling factor applied to the signal to test the sensitivity level of the 
experiment. Truncated Gaussian priors are used for each of the nuisance parameters Ok, which define the sensitivity of 
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the predicted signal and background estimates to systematic uncertainties. These can take the form of uncertainties 
on overall rates, as well as the shapes of the distributions used for combination. These systematic uncertainties can 
be far larger than the expected Higgs signal, and are therefore important in the calculation of limits. The truncation 
is applied so that no prediction of any signal or background in any bin is negative. The posterior density function is 
then integrated over all parameters (including correlations) except for R, and a 95% credibility level upper limit on 
R is estimated by calculating the value of R that corresponds to 95% of the area of the resulting distribution. 



B. Modified Frequentist Method 

The Modified Frequentist technique relies on the CL S method, using a log-likelihood ratio (LLR) as test statistic (22j: 

LLR=2ln PjMam (2) 

p(data|if ) 

where Hi denotes the test hypothesis, which admits the presence of SM backgrounds and a Higgs boson signal, while 
Hq is the null hypothesis, for only SM backgrounds. The probabilities p are computed using the best-fit values of the 
nuisance parameters for each event, separately for each of the two hypotheses, and include the Poisson probabilities 
of observing the data multiplied by Gaussian constraints for the values of the nuisance parameters. This technique 
extends the LEP procedure which does not involve a fit, in order to yield better sensitivity when expected signals are 



small and systematic uncertainties on backgrounds are large 24 1. 

The CL S technique involves computing two p-values, CL s +b and CLi,. The latter is defined by 

l-CL b = p(LLR < LLR ohs \H ), (3) 

where LLR \> S is the value of the test statistic computed for the data. 1 — CL& is the probability of observing a 
signal-plus-background-like outcome without the presence of signal, i.e. the probability that an upward fluctuation of 
the background provides a signal-plus-background-like response as observed in data. The other p-value is defined by 

CL s+h =p{LLR>LLR ohs \H x ) : (4) 

and this corresponds to the probability of a downward fluctuation of the sum of signal and background in the data. A 
small value of CL s+ b reflects inconsistency with H\. It is also possible to have a downward fluctuation in data even in 
the absence of any signal, and a small value of CL s+ b is possible even if the expected signal is so small that it cannot be 
tested with the experiment. To eliminate the possibility of excluding a signal to which there is insufficient sensitivity 
(an outcome expected 5% of the time at the 95% C.L., for full coverage), we use the quantity CL S — CL s+ b/CLb. If 
CL S < 0.05 for a particular choice of Hi, that hypothesis is deemed excluded at the 95% C.L. 

Systematic uncertainties are included by fluctuating the predictions for signal and background rates in each bin of 
each histogram in a correlated way when computing CL s +e, and CLb. 



C. Systematic Uncertainties 

The uncertainty on the measurement of the integrated luminosity is 5.8% (CDF) and 6.1% (DO). Of this value, 4% 
arises from the uncertainty on the inelastic pp scattering cross section, which is correlated between CDF and DO. The 
uncertainty on the rates for ti production and for single and di-electroweak boson production are taken as correlated 
between the two experiments. As the methods of measuring the multi-jet ("QCD") backgrounds differ between CDF 
and DO, there is no correlation assumed between these rates. The calibrations of fake leptons, unvetoed 7 — > e + e~ 
conversions, 6-tag efficiencies and mistag rates are performed by each collaboration using independent data samples 
and methods, hence are considered uncorrelated. 

Tables [TT] to IVIIII summarize the various contributions to the systematics uncertainties to the input distributions 
used in the limit setting, broken down by experiment and channel. Entries in the tables labeled as "shape" systematics 
do not have the same value across all bins of the relevant distribution and model the systematic variation of the shape 
for that source of uncertainty. In these cases the number given is the event weighted mean fluctuation away from the 
nominal distribution - i.e. related to the flat component of the uncertainty. 
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Contribution 


Signal Z 


— > e + e 


Z — > TT 


tt 


diboson 


QCD 


Jet energy scale (shape) 


0.12 


+0.30 
+0.22 


+0.05 

0.00 


0.56 


0.73 


0.0 


Electron identification 


2.4 


2.4 


2.4 


2.4 


2.4 


0.0 


Electron energy scale (shape) 


+0.32 

-0.23 


0.00 
+0.30 


+0.77 
-0.50 


— 0.19 
+0.28 


— 0.09 
+0.16 


0.0 


Tau identification 


4.2 


4.2 


4.2 


4.2 


4.2 


0.0 


Tau energy scale (shape) 


+0.29 
-0.22 


+0.22 
+0.23 


+0.82 
-0.63 


+0.42 
-0.75 


+0.63 
-0.31 


0.0 


gg — >• A acceptance 


2.1 


0.0 


0.0 


0.0 


0.0 


0.0 


bb — > A acceptance 


3.6 


0.0 


0.0 


0.0 


0.0 


0.0 


MC Cross sections 


0.0 


2.2 


2.2 


10.0 


6.0 


0.0 


QCD 


0.0 


0.0 


0.0 


0.0 


0.0 


15.0 


Luminosity 


5.8 


0.0 


5.8 


5.8 


5.8 


0.0 



TABLE II: Percentage systematic uncertainties for each distribution in the CDF r^rhad analysis. Signal uncertainties are for 
M A = 130 GeV. 



Contribution 


Signal Z 




Z -^TT 


tt 


diboson QCD 


Jet energy scale (shape) 
Muon identification 


0.07 
2.7 


+0.24 
-0.38 

2.7 


0.0 
2.7 


+0.54 
-0.48 

2.7 


+ U.46 
-0.58 
2.7 


0.0 
0.0 


Tau identification 


4.2 


4.2 


4.2 


4.2 


4.2 


0.0 


Tau energy scale (shape) 


+0.23 
-0.08 


0.0 


+0.54 
-0.77 


+0.97 
-0.75 


+0.40 
-0.70 


0.0 


gg — > A acceptance 


2.1 


0.0 


0.0 


0.0 


0.0 


0.0 


bb — > A acceptance 


3.6 


0.0 


0.0 


0.0 


0.0 


0.0 


MC cross sections 


0.0 


2.2 


2.2 


10.0 


6.0 


0.0 


QCD 


0.0 


0.0 


0.0 


0.0 


0.0 


20.0 


Luminosity 


5.8 


0.0 


5.8 


5.8 


5.8 


0.0 



TABLE III: Percentage systematic uncertainties for each distribution in the CDF r M Th a d analysis. Signal uncertainties are for 
M A = 130 GeV. 



Contribution 


Signal Z 


-+ e+e 


Z — S- 


Z -^TT 


tt 


diboson QCD 


Electron energy scale (shape) 
Jet energy scale (shape) 
Electron identification 


+0.23 

-0.46 
-0.08 
0.00 

2.4 


0.0 
0.0 
2.4 


-0.36 
-0.60 
-0.34 
-0.30 

2.4 


+0.72 
-0.62 
-0.05 
0.00 

2.4 


0.00 
+0.26 

0.57 
2.4 


0.0 
0.29 
2.4 


0.0 
0.0 
0.0 


Muon identification 


2.7 


2.7 


2.7 


2.7 


2.7 


2.7 


0.0 


gg — > A acceptance 


2.1 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


bb — > A acceptance 


3.6 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


MC Cross sections 


0.0 


2.2 


2.2 


2.2 


10.0 


6.0 


0.0 


QCD 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


20.0 


Luminosity 


5.8 


5.8 


5.8 


5.8 


5.8 


5.8 


0.0 



TABLE IV: Percentage systematic uncertainties for each distribution in the CDF t^t^ analysis. Signal uncertainties are for 
M A = 130 GeV. 

IV. COMBINED RESULTS 

Tables iLXl and IXl give the 95% confidence limits on the cross section x branching ratio for MSSM Higgs production 
and decay in the di-tau channel, using the two different approaches outlined above. Good agreement in the results 
for the two procedures is seen with variations at less than 10%. The results are shown graphically in Figure [3l using 
the CL S calculations. The observed limits are generally in good agreement with expectation with no evidence for 
significant excess for 90 < Ma < 200 GeV. 
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Contribution 


Signal diboson QCD 


tt w- 


-> Iv 


Z -> e+e~ 


Z -^tt 


Electron Identification 


3 


3 


3.3 


0.0 


3.3 


3 


3 


3.3 


3.3 


Electron-tau fake rate 








0.0 


0.0 


0.0 








13 


0.0 


Tau identification 


6 





5.3 


0.0 


7.1 


5 


6 


3.9 


4.1 


Tau track reconstruction 


1 





1.0 


0.0 


1.0 


1 





1.0 


1.0 


Tau energy scale (shape) 





4 


0.0 


0.0 


0.0 








0.0 


1.3 


Trigger (shape) 


3 


8 


4.1 


0.0 


3.0 


4 


4 


4.2 


5.9 


Signal acceptance 


4 





0.0 


0.0 


0.0 








0.0 


0.0 


MC cross sections 








5.0 


0.0 


5.0 








5.0 


5.0 


W+jets 








0.0 


0.0 


0.0 


6 


8 


0.0 


0.0 


QCD 








0.0 


13.0 


0.0 








0.0 


0.0 


Luminosity 


6 


1 


6.1 


0.0 


6.1 


6 


1 


6.1 


6.1 



TABLE V: Percentage systematic uncertainties for each distribution in the D0 r e Thad analysis - combined across all three tau 
categories. Signal uncertainties are for Ma = 130 GeV. 



Contribution 


Signal diboson QCD 


tt 


W -> iv Z - 




Z -^TT 


Muon identification 


1.1 


1.1 


0.0 


1.1 


1.1 


1.1 


1.1 


Tau identification 


4.2 


3.9 


0.0 


4.2 


5.6 


3.9 


3.9 


Tau track reconstruction 


1.0 


1.0 


0.0 


1.0 


1.0 


1.0 


1.0 


Tau energy scale (shape) 


0.79 


0.0 


0.0 


0.0 


0.0 


0.0 


1.3 


Trigger 


3.0 


3.0 


0.0 


3.0 


3.0 


3.0 


3.0 


Signal acceptance 


4.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


MC cross sections 


0.0 


5.0 


0.0 


5.0 


0.0 


5.0 


5.0 


W+jets 


0.0 


0.0 


0.0 


0.0 


13.0 


0.0 


0.0 


QCD 


0.0 


0.0 


32 


0.0 


0.0 


0.0 


0.0 


Luminosity 


6.1 


6.1 


0.0 


6.1 


6.1 


6.1 


6.1 



TABLE VI: Percentage systematic uncertainties for each distribution in the D0 r^Thad - (Runlla) - combined across all three 
tau categories. Signal uncertainties are for Ma = 130 GeV. 



Contribution 


Signal diboson 


Z — > 


Z tt 


QCD 


Muon identification 


4.0 


4.0 


4.0 


4.0 


4.0 


0.0 


Muon track reconstruction 


2.0 


2.0 


2.0 


2.0 


2.0 


0.0 


Tau identification 


3.9 


0.0 


0.0 


3.8 


0.0 


0.0 


Tau track reconstruction 


1.4 


0.0 


0.0 


1.4 


0.0 


0.0 


Tau energy scale 


2.6 


2.4 


2.7 


2.5 


2.3 


0.0 


Trigger 


5.0 


5.0 


5.0 


5.0 


5.0 


0.0 


Signal acceptance 


4.60 


0.0 


0.0 


0.0 


0.0 


0.0 


MC cross sections 


0.0 


5.0 


5.0 


5.0 


5.0 


0.0 


QCD 


0.0 


0.0 


0.0 


0.0 


0.0 


22 


Luminosity 


6.1 


6.1 


6.1 


6.1 


6.1 


0.0 



TABLE VII: Percentage systematic uncertainties for each distribution in the D0 r M rh a d - (Runllb) - combined across all three 
tau categories. Signal uncertainties are for Ma = 130 GeV. 

V. INTERPRETATION WITHIN THE MSSM 

Though at leading order the Higgs sector of the MSSM can be described with just two parameters, with higher 
order corrections comes a dependence on other model parameters. To interpret the exclusion within the MSSM these 
parameters are fixed in four benchmark scenarios (25| . The four scenarios considered are defined in terms of: Msusy > 
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Contribution 


Signal 


QCD 


tt 


W -> iv 


diboson Z 




Z — > 


Z -^tt 


Jet energy scale 


2.0 


0.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


Electron identification 


2.0 


0.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


Muon identification 


0.4 


0.0 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


Vertex modelling 


2.0 


0.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


Trigger 


4.0 


0.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


Signal acceptance 


4.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


MC cross sections 


0.0 


0.0 


5.0 


5.0 


5.0 


5.0 


5.0 


5.0 


QCD 


0.0 


7.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


Luminosity 


6.1 


0.0 


6.1 


6.1 


6.1 


6.1 


6.1 


6.1 



TABLE VIII: Percentage systematic uncertainties each distribution in the D0 r e r M analysis. Signal uncertainties are for Ma 
= 130 GeV. 





Observed 




Expected Limits / pb 




Mass / GeV 


Limits/pb 


-2a 


-la 


median 


+ la 


+2a 


90 


12.6 


10.6 


14.6 


20.2 


28.7 


39.2 


100 


17.4 


6.89 


9.52 


13.2 


18.6 


25.1 


110 


9.89 


3.30 


4.45 


6.18 


8.63 


11.6 


120 


4.22 


1.77 


2.38 


3.40 


4.80 


6.69 


130 


1.94 


1.16 


1.55 


2.16 


2.92 


3.93 


140 


1.41 


0.805 


1.08 


1.55 


2.10 


2.80 


150 


1.00 


0.615 


0.795 


1.13 


1.56 


2.14 


160 


0.832 


0.495 


0.655 


0.895 


1.26 


1.71 


170 


0.771 


0.405 


0.525 


0.745 


1.07 


1.41 


180 


0.647 


0.325 


0.445 


0.615 


0.865 


1.65 


190 


0.628 


0.295 


0.385 


0.545 


0.765 


1.04 


200 


0.629 


0.255 


0.335 


0.475 


0.665 


0.915 


TABLE IX: Combined Cross section x branching 


ratio limits using Bayes method. 




Observed 




Expected Limits / pb 




Mass / GeV 


Limits/pb 


-2a 


-la 


median 


+ la 


+2a 


90 


13.6 


11.0 


14.5 


20.0 


28.6 


39.3 


100 


17.7 


6.55 


8.77 


12.2 


17.2 


23.4 


110 


9.68 


3.21 


4.32 


5.76 


8.20 


11.1 


120 


4.19 


1.74 


2.39 


3.26 


4.60 


6.22 


130 


2.03 


1.07 


1.50 


2.14 


2.94 


4.03 


140 


1.45 


0.837 


1.04 


1.47 


2.08 


2.81 


150 


1.02 


0.586 


0.780 


1.11 


1.54 


2.09 


160 


0.829 


0.454 


0.633 


0.884 


1.24 


1.70 


170 


0.807 


0.406 


0.529 


0.719 


1.01 


1.38 


180 


0.697 


0.315 


0.431 


0.595 


0.841 


1.14 


190 


0.681 


0.281 


0.357 


0.514 


0.724 


1.00 


200 


0.653 


0.261 


0.325 


0.452 


0.638 


0.867 



TABLE X: Combined cross section x branching ratio limits using CLs. 

the mass scale of squarks, ^, the Higgs sector bilinear coupling, M 2 , the gaugino mass term, A t , the trilinear coupling 
of the stop sector, At,, the trilinear coupling of the sbottom sector and rrig the gluino mass term. The maximal- mixing, 
m™ 1 , scenario is defined as: 
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M A [GeV] 



FIG. 3: 95% Confidence limits on cross section x branching ratio. The solid black and dashed grey lines show the observed 
and expected limits respectively. The yellow and blue hatched bands around the expected limit show the 1 and 2cr deviations 
from the expectation. 



Msvsy = ITeV, [i = 200 GeV, M 2 = 200 GeV, 
X t = 2Msusy 
A b = A t , rrig = O.SMstjsy- 

and the no-mixing scenario - with vanishing mixing in the stop sector and a higher SUSY mass scale to avoid the 
LEP Higgs bounds: 



Msusy = 2TeV, fx = 200 GeV, M 2 = 200 GeV, 
X t = 0, A b = A t , m g = 0.8A/ SUSY . 

Four scenarios are constructed from these two by the consideration of both + and - signs for ji. 
Tables IXlH IXI| |XIV| and IXIIll give the observed and median expected 95% confidence limits on tan/3 for the tested 
mass hypotheses for the four different benchmark scenarios considered. This is shown graphically in Figure 2] 

In this preliminary result the signal cross sections and branching fractions within each scenario have been calculated 
using feynhiggs [4] - with gg — > H production from (26l - [36| and SM bb — >• H production from [33] and references 
therein and MRST2002 NNLO PDFs [38] - with no theoretical uncertainties considered. Tan/3 dependent width effects 
have not been included, though in the region of the tan/3- Ma plane where limits have been set these are not expected 
to strongly affect the limit Q. 
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Observed Expected Limits / pb 



M A GeV 


Limits 


-2a 


-la 


median 


+ la 


+2a 


90 


30 


27 


31 


36 


44 


51 


100 


44 


27 


31 


37 


44 


51 


110 


42 


24 


28 


32 


38 


44 


120 


34 


22 


25 


30 


35 


41 


130 


29 


21 


25 


30 


35 


40 


140 


29 


22 


25 


29 


35 


41 


150 


30 


23 


26 


31 


37 


43 


160 


32 


24 


28 


33 


39 


46 


170 


37 


27 


30 


35 


42 


49 


180 


41 


27 


32 


38 


45 


52 


190 


47 


30 


34 


41 


48 


56 


200 


53 


34 


38 


44 


52 


61 



TABLE XI: Combined 95% confidence limits on tan/3 for each mass hypothesis in the m h max and negative /i scenario. 



Observed Expected Limits / pb 



M A GeV 


Limits 


-2a 


-la 


median +la 


+2a 


90 


31 


28 


32 


37 


45 


53 


100 


46 


28 


32 


38 


45 


53 


110 


43 


25 


28 


33 


40 


46 


120 


34 


22 


26 


30 


36 


42 


130 


29 


21 


25 


30 


36 


42 


140 


30 


22 


25 


30 


36 


42 


150 


31 


23 


27 


32 


38 


44 


160 


33 


24 


29 


34 


40 


47 


170 


38 


27 


31 


36 


43 


50 


180 


42 


28 


33 


39 


46 


54 


190 


48 


31 


35 


42 


50 


59 


200 


55 


35 


39 


46 


54 


64 



TABLE XII: Combined 95% confidence limits on tan/3 for each mass hypothesis in the mh max and positive /j, scenario. 



This combination of Tevatron results from CDF and DO in the A — > rr channel sets the most stringent limits to 
date on the search for MSSM Higgs in that final state. 
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Observed Expected Limits / pb 



M A GeV 


Limits 


-2a 


-la 


median 


+ la 


+2a 


90 


30 


27 


31 


37 


44 


52 


100 


45 


27 


32 


37 


44 


52 


110 


42 


24 


28 


32 


38 


45 


120 


34 


22 


26 


30 


36 


41 


130 


29 


20 


25 


30 


35 


41 


140 


30 


23 


26 


30 


36 


42 


150 


30 


23 


26 


32 


37 


43 


160 


32 


24 


28 


33 


40 


46 


170 


38 


27 


31 


36 


42 


49 


180 


41 


28 


32 


38 


45 


52 


190 


47 


30 


34 


41 


49 


57 


200 


54 


34 


38 


45 


53 


62 



TABLE XIII: Combined 95% confidence limits on tan/3 for each mass hypothesis in the no-mixing and negative /i scenario. 



Observed Expected Limits / pb 



M A GeV 


Limits 


-2a 


-la 


median +la 


+2a 


90 


31 


27 


31 


37 


44 


52 


100 


45 


27 


32 


37 


44 


52 


110 


42 


24 


28 


32 


39 


45 


120 


34 


22 


26 


30 


36 


42 


130 


29 


20 


25 


30 


35 


41 


140 


30 


23 


26 


30 


36 


42 


150 


30 


23 


27 


32 


37 


43 


160 


33 


24 


28 


34 


40 


47 


170 


38 


27 


31 


36 


42 


50 


180 


41 


28 


32 


38 


45 


53 


190 


47 


31 


35 


41 


49 


58 


200 


54 


34 


38 


45 


53 


62 



TABLE XIV: Combined 95% confidence limits 



on tan/3 for each mass hypothesis in the no-mixing and positive /i scenario. 
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O2.100F 

c 

- 90 
80 
70 
60 
50 
40 
30 
20 
10 



Tevatron Run II Preliminary, L= 1.8-2.2 fb" 1 
m h max, n=-200 GeV 



I I Excluded by LEP 
— Observed limit 
Expected limit 
Expected limit ± 1 o 
l^^l Expected limit ± 2 o 




180 200 

m A [GeV/c 2 ] 



c^100 
c 

- 90 
80 
70 
60 
50 
40 
30 
20 
10 



Tevatron Run II Preliminary, L= 1.8-2.2 fb" 1 

m h max, n=+200 GeV 



Excluded by LEP 
- Observed limit 
Expected limit 
Expected limit lie 
] Expected limit±2o 




m A [GeV/C 1 ] 



cq.100f 
c 

*" 90: 
80 
70 
60 
50 
40 
30 
20 
10 



Tevatron Run II Preliminary, L= 1.8-2.2 fb" 1 
no mixing, u=-200 GeV 



Excluded by LEP 
■ Observed limit 
Expected limit 
Expected limit ± 1 a 
3 Expected limit ± 2 a 




180 200 

m A [GeV/c 2 ] 



ci100 
c 

- 90 
80 
70 
60 
50 
40 
30 
20 
10 



Tevatron Run II Preliminary, L= 1.8-2.2 fb" 1 
no mixing, u=+200 GeV 



I I Excluded by LEP 
— Observed limit 
Expected limit 
Expected limit ± 1 a 
t\^t Expected limit ± 2 a 




180 200 

m A [GeV/c 2 ] 



FIG. 4: 95% Confidence limits in the tan/3-MA plane for the 4 benchmark scenarios: maximal mixing (top) and no mixing 
(bottom) for fi < (left) and /i > (right). The black line denotes the observed limit, the grey line the expected limit and the 
hatched yellow and blue regions denote the ±1 and 2 a bands around the expectation. The shaded light-green area shows the 
limits from LEP. 
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